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Glass is everywhere. We use and are surrounded by glass objects which make tangible the reality 
of glass as a distinct state of matter. Yet, glass as we know it is usually obtained by cooling a liquid 
sufficiently rapidly below its melting point to avoid crystallisation. The viscosity of this supercooled 
liquid increases by many orders of magnitude upon cooling, until the liquid becomes essentially 
arrested on experimental timescales below the “glass transition” temperature. From a structural 
viewpoint, the obtained glass still very much resembles the disordered liquid, but from a mechanical 
viewpoint, it is as rigid as an ordered crystal. Does glass qualify as a separate state of matter? We 
provide a pedagogical perspective on this question using basic statistical mechanical concepts. We 
recall the definitions of states of matter and of phase transitions between them. We review recent 
theoretical results suggesting why and how an “ideal glass” can indeed be defined as a separate 
equilibrium state of matter. We discuss recent success of computer simulations trying to analyse 
this glass state. We close with some experimental perspectives. 


I. WHAT IS A GLASS? 


We possess an intuitive notion of what a state of mat- 
ter is. For a pure substance, there are usually three 
such states: gas, liquid, and crystal [1]. The gas is 
found at high temperatures and transforms into a liq- 
uid upon cooling. For an infinitely slow cooling process, 
where thermal equilibrium is maintained all along, con- 
densation at the boiling point is a first-order equilibrium 
phase transition associated with a discontinuous jump in 
various thermodynamic quantities, such as the specific 
volume shown in Fig. 1 (blue curve). The liquid itself 
transforms abruptly into a crystalline solid when cooled 
below the melting temperature T,, through another first- 
order equilibrium phase transition. These states of mat- 
ter are characterised and distinguished by their macro- 
scopic properties. For example, the gas occupies the en- 
tire available volume, unlike the crystalline solid which 
retains its shape. 

Glass is absent from this description because we have 
so far focused on equilibrium properties. If a liquid is 
cooled quickly enough below T,, (red curve in Fig. 1), 
crystallisation can be avoided and the system becomes a 
supercooled liquid [2]. This supercooled state can exist 
over very long times because the formation of a crystal is 
very slow for a large number of substances. In this super- 
cooled regime, the system appears for all practical pur- 
poses at thermal equilibrium although it is metastable 
with respect to the crystal. For instance, its physical 
properties do not depend on time and can easily be re- 
produced in independent experiments. 

While the structure of the supercooled liquid changes 
very little with decreasing the temperature, its trans- 
port properties, such as viscosity or diffusivity, change 
by many orders of magnitude over a small temperature 
range. At the microscopic scale, molecular motion slows 
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FIG. 1. Phases of matter. A gas transforms into a liq- 
uid below the boiling temperature Tj, through an equilibrium 
first-order phase transition. A liquid crystallises below the 
melting temperature Tm through another first-order phase 
transition. A metastable supercooled liquid can persist be- 
low Tm and gradually transforms into a solid at the glass 
transition temperature T,. Below T,, the glass is frozen in a 
rigid disordered structure which depends on the cooling rate. 
This is the common glass, which slowly ages towards equilib- 
rium. An equilibrium “ideal glass” state would be obtained 
across a first-order phase transition at Tk provided thermal 
equilibrium could be maintained down to Tx. 


down dramatically and becomes so slow that very little 
motion can be detected over timescales of several hun- 
dreds of seconds when the temperature falls below the 
glass transition temperature T,. Below T,, the super- 
cooled liquid is out of equilibrium, and molecules stop 
diffusing over experimentally accessible timescales [3]. 


This nearly-arrested liquid behaves macroscopically like 
a solid even though its microscopic structure still resem- 
bles the one of a completely disordered liquid. The “glass 
transition” temperature at T, is not a well-defined tem- 
perature and depends on the protocol used to cool the 
liquid. As a consequence, the glass that we use in our 
everyday life is a non-equilibrium material with physi- 
cal properties that depend on its preparation history, see 
Fig. 1, and may vary with time as the system slowly ages 
towards equilibrium. Therefore, a conventional glass is 
not a state of matter, as it is not an equilibrium mate- 
rial with crisply defined and fully reproducible thermo- 
dynamic properties. 

This is of course not the end of our story. The idea that 
a glass state of matter can exist in equilibrium conditions 
has a long history, dating back at least to the work of 
Kauzmann in 1948 [4]. Kauzmann collected calorimet- 
ric data measured at equilibrium above 7, for a num- 
ber of molecular liquids, and by extrapolation below T, 
noted the possibility that a thermodynamic phase tran- 
sition could take place at a temperature Tk < T,, now 
called the Kauzmann temperature. Although Kauzmann 
himself favoured an interpretation in terms of a “pseudo- 
critical point” marking the limit of stability of the liquid 
with respect to crystallization, his observation raised the 
prospect that the system below Tx could be in an equi- 
librium glass state. We shall call this state an “ideal 
glass”, to make the distinction with conventional non- 
equilibrium glasses very clear. The Kauzmann tempera- 
ture and the equilibrium ideal glass phase are shown in 
Fig. 1. 

The hypothesis that an equilibrium phase transition 
at Tx underlies the formation of conventional glasses has 
given rise to a large litterature [3]. Our understanding of 
its nature has significantly evolved over the years since 
the early work of Gibbs and diMarzio who first described 
a simple lattice model with a thermodynamic liquid-glass 
transition [5]. Here, we offer a pedagogical discussion of 
the ideal glass and the equilibrium glass transition relying 
on recent theoretical and numerical developments. Our 
aim is to make this discussion accessible to anyone with 
an undergraduate background in statistical mechanics at 
the level of well-known textbooks [6-8]. 

In Sec. II, we review the concepts of states of matter 
and phase transitions to precisely define the question ad- 
dressed in our title. In Sec. II, we summarise the mean- 
field theory of the glass transition and explain why, in 
this limit, glass does correspond to a state of matter. 
In Sec. IV, we present numerical results of glass-forming 
liquids which are consistent with the existence of an un- 
derlying liquid-to-glass equilibrium phase transition. We 
conclude in Sec. V. 


II. WHAT IS A STATE OF MATTER? 


Following the conventional Landau approach to phase 
transitions [6] and without delving into all the subtleties, 


a state or a phase of matter is commonly defined as a 
homogeneous region of space in which all intensive ther- 
modynamic variables are well-defined and uniform. We 
distinguish among these variables the control parameters 
(temperature, pressure, magnetic field, etc.) which can 
be imposed by an external operator, and others, which 
are densities associated with extensive quantities (spe- 
cific mass, magnetisation, massic energy, etc.). The or- 
der parameter is one of these variables, but it plays a 
special role since it takes different values in the different 
phases, and its behaviour makes it possible to distinguish 
between them. 


The different phases are usually represented by do- 
mains in a phase diagram, which is the space of control 
parameters. These domains are limited by phase tran- 
sition lines through which the order parameter can vary 
continuously (second-order phase transitions) or discon- 
tinuously (first-order phase transitions). To determine 
the domain of stability of the different phases, Landau 
proposed a systematic method based on statistical me- 
chanics [6], which now corresponds to the starting point 
of most introductory courses on phase transitions [7]. 
Having identified a good order parameter, one has to 
express the equilibrium free energy of the system as a 
function of this order parameter for fixed values of the 
control parameters (for example the temperature T), and 
then determine its absolute minimum for each value of 
T. The computed values of the order parameter and of 
the free energy at this minimum represent the value that 
would be measured in a macroscopic system at thermal 
equilibrium. 

Some phase transitions are associated with a sponta- 
neous symmetry breaking, corresponding to a change in 
the degree of symmetry of the system [9]. The varia- 
tion of the order parameter from one phase to another 
reflects this loss of symmetry. For example, in the case 
of the transition between paramagnetic and ferromag- 
netic phases, the order parameter is the magnetisation 
vector which has a non-zero value in the ferromagnetic 
phase. The magnetisation then favours a single direction 
of space and thus breaks rotational invariance. Similarly, 
the crystallisation of a liquid into an ordered periodic 
lattice decreases the degree of symmetry of the system. 
Some phase transitions are not associated with a symme- 
try breaking, such as the liquid-gas transition. 


To describe glass as an equilibrium phase along similar 
lines, a legitimate objection emerges immediately. Glass 
as we know it refers to a liquid that has left thermal equi- 
librium at J, and can no longer explore its configuration 
space ergodically. This means that an equilibrium free 
energy cannot even be defined for a conventional glass 
and the above discussion appears irrelevant. Therefore, 
the questions we wish to address are: assuming thermal 
equilibrium can be maintained at arbitrary low tempera- 
tures and configuration space can be explored in equilib- 
rium conditions, does the free energy show a transition 
from a liquid to a glass state, what is the order parameter 
that would then reveal the corresponding phase transi- 


tion, and what is the corresponding symmetry breaking? 


III. LANDAU THEORY OF THE GLASS 
TRANSITION IN THE MEAN-FIELD LIMIT 


A series of works initiated in 2012 has led to the com- 
plete derivation of an analytic solution of the glass transi- 
tion for simple models of glass-forming liquids which be- 
comes mathematically exact in the limit of a large num- 
ber of space dimensions, d — oo [10]. This non-physical 
mathematical limit makes the statistical mechanics of su- 
percooled liquids exactly solvable. Such a situation is 
very common in statistical mechanics where it is often 
found that mean-field approximate solutions in finite di- 
mensions become exact when d becomes large [7]. In 
large dimensions, the interactions between a given par- 
ticle and its neighbours can be replaced by an average 
force which can be self-consistently determined. De- 
spite this immense simplification, the calculation remains 
mathematically involved because it requires the explicit 
description of the aperiodic structure characterising the 
glass phase. The results found in large dimensions con- 
firm a long series of earlier results which hinted at the 
nature of the glass transition when using mean-field ap- 
proximations [11—13). 

It is possible to outline the main results obtained in 
this framework following the Landau approach described 
in Sec. II. To this end, the first step is to identify the good 
order parameter, Q, which can distinguish between the 
equilibrium liquid and glass phases, and then to express 
the free energy of the system as a function of Q, V(Q). 
In this context, a suitable choice for Q is what is called 
the overlap and for the associated V(Q) the Franz-Parisi 
potential [14]. 

What indeed is a good order parameter Q? Clearly, 
neither the average density nor its spatial fluctuations 
can be used because if the glass looks like an arrested 
liquid then its density, p, its pair correlation function 
g(r) and its static structure factor S'(k) also look similar. 
Therefore, contrary to the description of the liquid state 
where they play a key role [15], density fluctuations are 
uninteresting spectators of the glass transition. 

The correct answer originates from another physics 
problem where glassy behaviour is also found, namely, 
spin glasses [16, 17]. Although spin glass physics has 
a much shorter history than structural glass physics, 
the formulation of statistical mechanics models for spin 
glasses and a deep analytic understanding of their prop- 
erties went comparatively much faster. Spin glass mod- 
els describe spins with interactions mediated by random 
coupling constants, which makes it possible to describe 
the properties of dilute magnetic alloys [18]. The appro- 
riate order parameter to distinguish the spin glass and 
paramagnetic phases is the overlap Q between magneti- 
sation profiles of two independent copies (usually called 
“replicas” in the field of disordered systems) of the same 
system [18]. For Ising spins, this is simply expressed as 
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FIG. 2. Overlap order parameter Q of the liquid-glass 
transition. It quantifies the degree of similarity between 
the density profiles of two independent equilibrium configu- 
rations (here obtained from numerical simulations of a two- 
dimensional glass-forming liquid [19]). For T > Tx, the liquid 
explores a large number of equilibrium structures and Q ~ 0. 
Instead, in the equilibrium glass for T < Tx, all equilibrium 
configurations are very similar and Q > 0. 


Q = (1/N) ON, 85, where ${) = +1 is the value 


a 


of spin i = 1,...,N in one copy and 3?) the value of 
the same spin in the second copy. Clearly, the overlap 
Q ~ 1 if copies 1 and 2 are highly correlated (spin glass), 
whereas for uncorrelated configurations Q ~ 0 (paramag- 
net). The overlap Q then quantifies the degree of similar- 
ity of the local degrees of freedom (here, the spins) in two 
configurations of the system drawn from the equilibrium 
Boltzmann distribution. 


For supercooled liquids, the overlap quantifies the de- 
gree of similarity between the density profiles measured 
in two independent configurations. As sketched in Fig. 2, 
for two copies in which the density profiles are uncorre- 
lated, Q ~ 0 but the overlap becomes close to 1 for cor- 
related configurations. The order parameter Q can also 
be interpreted as a distance between the two copies in 


configuration space [20]. 

In practice the calculation of Q requires the introduc- 
tion of a microscopic length scale a corresponding to the 
amplitude of the thermal vibrations around the average 
positions of the particles, roughly corresponding to the 
cage size [21]. If one considers two copies of the liquid 
Ca = GPO ai ay defined from the particle positions 


ri) for a = 1, 2, the microscopic density profiles are de- 


fined as p( (r) = Al 6(r — rv‘) with 6(x) the Dirac 
p i=1 
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distribution. The overlap then reads 


QlC Ca) = = ff arar’p (r)p(r"yw( lr =r") 
1 & (1) 


with w(x) = 1 if x < a and 0 otherwise. From Eq. (1) 
it is clear that Q[C1,C2] ~ 0 for uncorrelated configura- 
tions, while Q[C1,C2] ~ 1 for strongly correlated ones. 
Just as for spin glasses, it is an intensive, normalised 
function which quantifies the degree of similarity of the 
main microscopic degrees of freedom in two independent 
configurations. 

The calculation of the Landau free energy V(Q) pro- 
ceeds in two separate steps. First, one chooses and fixes a 
reference configuration Co for the system at thermal equi- 
librium. The overlap between Cg and any other configu- 
ration C is Q|C,Co], defined by Eq. (1). The equilibrium 
probability distribution function P(Q,Co) of this overlap 
is obtained by performing a canonical average 


P(Q,Co) = > Pea(C)5 (Q — QIC,Co]), (2) 
Cc 


where the statistical weight of any configuration C is 
controlled by the Boltzmann factor at temperature T 
(kpT = 1/8), 


e~BHIC] 
Peq(C) = Yio, BHC (3) 


By definition, the free energy V(Q,Co) is the thermody- 
namic potential controlling the fluctuations of the overlap 
Q for a fixed reference configuration [6], 


P(Q,Co) ~ exp(—BNV(Q,Co)), (4) 


so that the most probable value of the overlap indeed 
corresponds to the minimum of the free energy in the 
thermodynamic limit, N — oo. In a second step one 
needs to perform a final average over all possible choices 
for the reference configuration to get the desired averaged 
free energy 


V(Q) = >- Pea(Co)V (Q, Co), (5) 


Co 


where the statistical weight of a given reference con- 
figuration is again controlled by the Boltzmann factor. 


The amorphous structure of the glass directly enters the 
canonical average in Eq. (2) which contains quenched dis- 
order (the reference configuration Cg) requiring a sepa- 
rate average in Eq. (5). It is this double averaging pro- 
cedure which makes the analytic calculation of V(Q) a 
difficult mathematical problem, even in the limit of infi- 
nite dimensions [10]. 


The temperature evolution of V(Q) in the mean-field 
limit is shown in Fig. 3(a). At sufficiently high tempera- 
tures, the free energy has a unique minimum for Q = 0. 
This equilibrium state corresponds to the liquid phase in 
which there is a large number of distinct accessible con- 
figurations, so that two randomly chosen configurations 
are very different. At lower temperature, the free energy 
develops a secondary minimum for Q@ > 0, which there- 
fore represents a metastable phase. This corresponds to a 
phase with strong correlations between equilibrium con- 
figurations which are all amorphous. When the tempera- 
ture decreases further, the secondary minimum becomes 
more and more stable. When the Kauzmann tempera- 
ture Tx is reached, the two minima have the same free 
energy and an equilibrium phase transition takes place 
between the liquid with Q = 0 and the ideal glass at 
Q > 0. This is a first-order phase transition because 
the overlap Q jumps discontinuously at Tx, as shown in 
Fig. 3(c). The central conclusion for our discussion is 
that, in the mean-field limit, an equilibrium glass exists 
as a genuine state of matter for T < Tx. 


The discontinuity of the order parameter Q at the tran- 
sition indicates that equilibrium configurations in the 
ideal glass phase are strongly correlated. More funda- 
mentally, the transition reflects a spontaneous symmetry 
breaking, which is that of the underlying replica sym- 
metry, and the discontinuity is associated with a spe- 
cific breaking scheme of this symmetry [16]. The liquid 
is symmetric with respect to the replica symmetry, be- 
cause if several replicas (copies) are sampled in the liquid, 
no density field is privileged. Conversely, the ideal glass 
breaks this symmetry, as now all replicas possess a similar 
amorphous density profile. Physically, the ideal glass is a 
disordered state of matter in real space (the distribution 
of the average positions of the particles is amorphous), 
but ordered in configuration space (all equilibrium con- 
figurations look similar). 


For T > Tx, the free energy difference AV(T) between 
the two minima of V(Q) represents the thermodynamic 
cost to force the system to remain localised near a given 
configuration in order to maintain a high value of Q, pre- 
venting the system to explore the much larger configura- 
tion space available in the liquid. The free energy cost 
is therefore entropic in nature. This free energy differ- 
ence AV(T) in fact represents a natural definition of the 
configurational entropy Scon¢(T) = AV(T)/T with the 
consequence that it vanishes continuously when T > Tx 
as shown in Fig. 3(b). This last result implies that the 
liquid-glass transition has no latent heat as the latent 
heat is normally proportional to the entropy jump at a 
first-order phase transition. The ideal glass transition is 
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FIG. 3. Mean-field theory of glass formation. (a) The Landau free energy V(Q) develops a secondary minimum corre- 
sponding to the metastable glass phase approaching Tx from above. (b) The free energy difference between the stable liquid 
and the metastable glass defines the configurational entropy Sconr(T'), which vanishes continuously at Tx, (c) where the overlap 
order parameter Q jumps discontinuously to a finite value. (d) Increasing the thermodynamic field € conjugate to Q induces a 
first-order transition line for T > Tx ending at a critical point in the universality class of the random-field Ising model. Above 
the critical point, the line of maximum overlap fluctuations is marked as dash-dotted line. 


thus an unconventional example of a discontinuous tran- 
sition [22]. More strikingly perhaps, the continuous van- 
ishing of the configurational entropy in Fig. 3(b) exactly 
realises the scenario in which the ideal glass transition 
corresponds to the entropy crisis first put forward by 
Kauzmann [4]. With hindsight, Kauzmann’s work ap- 
pears truly visionary. 


This exposition of the mean-field theory of the glass 
transition closely follows traditional descriptions of the 
van der Waals model for the liquid-gas transition or the 
Curie-Weiss model to describe ferromagnetism. Push- 
ing these analogies further, one can introduce the field 
€ that is thermodynamically coupled to the order pa- 
rameter Q [14]. This field plays a role similar to the 
magnetic field H conjugate to the magnetisation M in 
the Curie-Weiss model. It introduces a new contribution 
—eQ in the Hamiltonian which for « > 0 makes large Q 
values energetically more favourable. In the presence of a 
field, equilibrium is found by looking for the minimum of 
the tilted free energy V(Q) — eQ for fixed (T, 6). When 


the glass is metastable for T > Tx, increasing € trig- 
gers a first-order phase transition towards a high-overlap 
state, see Fig. 3(d). Physically, this can be understood 
from an energy-entropy argument, as for many equilib- 
rium phase transitions [9]. When « is increased, the sys- 
tem gains energy by adopting a large value of Q but 
loses configurational entropy because it only visits a set 
of highly-correlated configurations. Beyond a threshold 
value €*(T) ~ TScont(T) = AV(T), the system ends up 
in a high-overlap phase. As a result, the existence of a 
discontinuous phase transition at the Kauzmann temper- 
ature Tx implies the existence of a first-order transition 
line in the (7, €) plane. This line ends at a second-order 
critical point, once again in full analogy with the liquid- 
gas transition in the van der Waals description. 


The behaviour of the order parameter in the vicin- 
ity of a critical point defines universality classes. Each 
universality class is characterised by a different set of 
critical exponents. For glasses, the presence of disor- 
der in the various canonical averages in Eqs. (2, 5) im- 


plies that the critical end-point at a finite applied field 
e belongs to the universality class of the random-field 
Ising model [23, 24]. (See Ref. [25] for an introduction to 
the random-field Ising model.) By contrast, the critical 
points of the liquid-gas and paramagnetic-ferromagnetic 
phase transitions both belong to the universality class 
of the Ising model with no disorder. The transition be- 
tween a paramagnet and a spin glass is also controlled 
by the presence of quenched disorder but it forms a uni- 
versality class on its own [16]. Despite the mathematical 
similarities between molecular glasses and conventional 
spin glasses, in particular the definition of the order pa- 
rameter, the emergence of a glassy phase is qualitatively 
distinct in these two systems. 

Before closing this section, it is instructive to pause and 
ask: should one be satisfied that after decades of theoret- 
ical work on the problem of the glass transition, the only 
solid set of results is about the mean-field theory? There 
are several answers. First, the study of phase transi- 
tions and critical phenomena in the presence of quenched 
disorder is notoriously more complicated than for pure, 
conventional systems [26]. It should thus come as no 
surprise that the statistical mechanics of glasses is a very 
difficult analytic problem, and the existence of a solid 
starting point is an important step. However, past ex- 
perience with disordered systems also suggests that in- 
cluding corrections to mean-field is going to be an even 
harder problem. At present it is not known whether and 
how the mean-field construction can survive the inclu- 
sion of fluctuations that are neglected in the mean-field 
theory. It is conceivable that the very existence of an 
ideal glass as a phase of matter, established within the 
mean-field approximation, becomes impossible when the 
relevant fluctuations are understood and included in the 
calculations. 


IV. DOES THE IDEAL GLASS EXIST BEYOND 
MEAN-FIELD? 


The ideal glass, which is defined in full equilibrium con- 
ditions, differs from the conventional material prepared 
in the laboratory that has fallen out of equilibrium near 
T, while cooling. How can one observe an ideal glass? 
The most direct way to detect the ideal glass consists in 
the quasi-static cooling of two copies of the supercooled 
liquid to maintain thermal equilibrium. A discontinuous 
jump of their relative overlap Q would be observed at the 
Kauzmann temperature Tx. This is conceptually very 
simple! This seems however impossible due to the dra- 
matic slowing down of the dynamics of the supercooled 
liquid at low temperature which prevents an ergodic ex- 
ploration of the configuration space below T,. 

The difference between ideal and conventional glasses 
appears subtle: the conventional glass is confined in con- 
figuration space because it moves too slowly to explore 
the many different available states, while the ideal glass is 
confined in configuration space because there is nowhere 


to go even given an infinite amount of time. Establish- 
ing this difference for temperatures where equilibration 
is extremely long is then obviously a challenge. In addi- 
tion, the first-order nature of the equilibrium liquid-glass 
transition implies that it cannot be simply revealed by 
the growth of some critical fluctuations, as would be the 
case when approaching a second-order phase transition. 
More subtle signs of an incipient phase transition need 
to be detected. 

Major progress has recently been made to tackle these 
two distinct problems in computer simulation studies. To 
solve the first problem, efficient Monte Carlo algorithms 
have been developed which are able to easily sample the 
configuration space in equilibrium conditions employing 
non-physical particle motions [27-31]. As a result, equi- 
librium thermodynamic measurements can be performed 
at low temperatures where this task would be impossible 
by simply simulating the physical dynamics. In essence, 
these algorithms solve the first problem and maintain 
equilibrium at temperatures much lower than in conven- 
tional simulations. For certain numerical models, it is es- 
timated that temperatures smaller than the experimental 
glass temperature J, can be equilibrated, thus provind- 
ing access to the experimental temperature gap between 
Te and Tr. 

The second problem concerns numerical signatures of 
the approach to a first-order phase transition. The 
most direct measurement amounts to directly estimating 
the Landau-Franz-Parisi thermodynamic potential V(Q) 
across a range of temperatures. This program was started 
long ago [32, 33] but was plagued by both equilibration 
problems and sampling issues. More recently, improved 
methods to assess V(Q) from Eq. (2) based on large de- 
viations techniques have been proposed [34] and com- 
bined with enhanced sampling techniques [35]. As a re- 
sult, the Landau free energy V(Q) has been measured 
in several numerical models of glass-forming liquids in 
dimension d = 2 and d = 3 down to very low tempera- 
tures [19, 34, 36, 37]. The expected behaviour for V(Q) 
has been confirmed directly in these computer simula- 
tions, as illustrated in Fig. 4(a), with the emergence of 
singularities at large Q when T is low enough. Since 
these measurements are performed in equilibrium condi- 
tions and do not require any extrapolation, they are di- 
rect proof that the system is on its way to an equilibrium 
first-order phase transition with an emerging metastable 
glass phase at large Q. 

Glass metastability can be established even more 
clearly by coupling the overlap Q to its conjugate field 
€ [19, 34, 36, 38, 39]. In three dimensions, one finds a 
discontinuous jump of the overlap as € is increased at 
constant T beyond a critical value «*(T'), see Fig. 4(c). 
This discontinuous jump disappears when the tempera- 
ture is too large, as expected from mean-field theory. 

As a result, a line of first-order transition «*(T') can be 
determined in the phase diagram (T, €), see Fig. 4(d), 
which ends at a critical point whose universality class 
was recently analysed by detailed finite-size scaling anal- 
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FIG. 4. Numerical results of three-dimensional systems reveal features of the mean-field theory. (a) Temperature 
evolution of Landau free energy V(Q) which develops a singularity at large Q consistent with the existence of a metastable 
glass phase. The dashed line highlights the non-convexity of the free energy at the lowest temperature. (b) Configurational 
entropy measured in two glass-forming models. The free energy difference between the equilibrium liquid and metastable glass 
decreases rapidly with decreasing T and seems to extrapolate to a finite Kauzmann temperature Tx. (c) The growth of the 
overlap upon application of a field € is smooth at large T’ but becomes discontinuous at low 7’, revealing a first-order equilibrium 
transition line e*(T). (d) Equilibrium phase diagram (T’, «). Dashed lines are temperatures shown in (a, c). The dash-dotted 
line represents the location of the inflection point of Q(e) which transforms in a discontinuity (full line) below a second-order 


critical point (dot). 


ysis [38, 40]. (At the critical point, the amplitude of 
the order parameter fluctuations only diverges in the 
thermodynamic limit, so that analysis of this divergence 
with system size yields information on the critical expo- 
nents [41].) Results in three-dimensional systems appear 
compatible with the predicted random-field Ising model 
universality class. 


Interestingly, the connection with the random-field 
Ising model implies that neither the critical point nor the 
first-order transition line can exist in two dimensions at 
any finite temperature. This is because the lower critical 
dimension [9] of the random-field Ising model is d = 2, 
so that for d < 2 there is no longer a finite transition 
temperature in this model [42]. This non-trivial predic- 
tion was recently confirmed in computer simulations of a 
two-dimensional glass-forming model [19]. 


Using similar techniques [43], the configurational en- 


tropy Scont(T’) was recently estimated in computer mod- 
els down to extremely low temperatures, including below 
the experimental T,, see Fig. 4(b). These measurements 
on the one hand confirm and amplify the results obtained 
from calorimetric experiments, such as the results dis- 
cussed in Kauzmann’s original work [4]. Whereas the 
nature of the configurational entropy obtained experi- 
mentally is subject to discussions [44], in the simulations 
Scont(T) can be interpreted as the free energy difference 
between the equilibrium liquid and metastable glass. The 
steep decrease observed in Fig. 4(b) at low temperatures 
shows that the glass could become the equilibrium phase 
below a finite temperature Tx. Two-dimensional simu- 
lations [45] are compatible with a value Tk = 0, again 
suggesting that d = 2 plays a special role with a zero- 
temperature critical point and no equilibrium glass phase 
at any T > 0. 


Finally, metastability near a first-order phase transi- 
tion is generically associated with an important length 
scale related to the physics of nucleation. If a system is 
prepared in a metastable phase, the emergence of the sta- 
ble phase proceeds by first nucleating a critical droplet of 
the stable phase which can then invade the entire system. 
The minimal size of the nucleating droplet is inversely 
proportional to the free energy difference between the two 
phases and diverges at the transition. This behaviour is 
generic to first-order phase transitions and describes for 
instance how crystals form below T,,. In the present case, 
one would essentially like to estimate the size of the criti- 
cal nucleus of the stable liquid from the metastable glass 
above Tx. However, considering the unconventional first- 
order nature of the transition at Tx, this is far from trivial 
and special computational tools are needed. Such tools 
have been proposed [46] and implemented [37, 47-49], al- 
lowing one to estimate this important length scale, €(T), 
which is also known as the point-to-set length scale [50]. 

Numerical measurements show that €(T') increases 
when the temperature 7 decreases, as expected when 
approaching a first-order phase transition. The tem- 
perature dependence of €(T) is compatible with the 
above nucleation argument, which suggests that €(T’) ~ 
1/Scont(T). The growth of the point-to-set length scale 
&(T) therefore accompanies the steep decrease in the con- 
figurational entropy which signals that glass and liquid 
free energies are getting closer as T decreases towards 
Tk. 


V. CONCLUSION 


State-of-the-art computer studies provide solid confir- 
mations that supercooled liquids carefully maintained in 
equilibrium conditions at low temperatures (even below 
T,) display thermodynamic fluctuations that are fully 
compatible with the mean-field description of an incip- 
ient discontinuous liquid-glass phase transition. These 
observations invalidate the well-known motto that that 
glass formation is essentially a dynamic process accom- 
panied by little structural change. Of course, the exis- 
tence of the transition at a finite Kauzmann temperature 
Tx > 0 remains to be demonstrated by direct measure- 
ments. This would represent a major objective for future 
work. From a computational viewpoint, this requires the 
development of numerical algorithms more sophisticated 
than current ones in order to achieve equilibration even 
closer to the putative Tk. 

What about experiments? Seventy years after Kauz- 
mann [4], calorimetric measurements provide estimates 
of a configurational entropy which robustly extrapolate 
to a finite Kauzmann temperature [44]. In recent years, 
indirect experimental estimates of the point-to-set corre- 
lation length scale €(T’) were also obtained via non-linear 
dielectric susceptibilities [51, 52]. (If the frequency of the 
imposed oscillatory electric field is larger than the typical 


relaxation frequency of the supercooled liquid, then the 
liquid approximately behaves like a mosaic of indepen- 
dent glasses of typical size €(T’) [22], each one of them re- 
sponding with an induced dipolar moment that depends 
on their typical size €(T), which can then be inferred.) 

The seemingly impossible task to achieve equilibration 
below T,, was overcome after the discovery that amor- 
phous glassy films prepared by physical vapour depo- 
sition are much closer to equilibrium than their liquid- 
cooled counterparts [53]. When the deposition rate is 
sufficiently low (of the order of 1 nm/s), and the tem- 
perature judiciously chosen (about 15 % lower than the 
experimental glass transition temperature T,), it has be- 
come possible to create glassy films that seem to remain 
in equilibrium much closer to Tx than in any previous 
work [54], thus providing novel support for previous ex- 
trapolations of the configurational entropy from above 
T,. These new glassy materials are “ultrastable”, be- 
cause they correspond to thermodynamic conditions im- 
possible to achieve using bulk cooling [55]. 

In addition to providing equilibrium materials much 
closer to Tx, the transformation kinetics of an ultrastable 
glass into a liquid upon heating closely resembles the nu- 
cleation and growth process commonly observed near a 
first-order phase transition [56, 57]. For thin ultrastable 
films, the transformation is initiated at the free surface 
and propagates as a front. This is reminiscent of hetero- 
geneous nucleation. For thicker samples, liquid droplets 
first nucleate in the bulk and their size then increases 
until complete transformation, reminiscent of homoge- 
neous nucleation [58]. These analogies reinforce the idea 
that ultrastable glasses behave very much like a distinct 
phase of matter accessed via a discontinuous phase tran- 
sition. Computer simulations of this transformation pro- 
cess [59, 60] could help unravel the main mechanisms at 
play, and could allow direct quantitive comparisons with 
predictions from classical nucleation theory [61]. 

The construction of a complete theoretical description 
of the nature of the glass remains an active research area. 
The progress summarised here results from a synergy be- 
tween mathematical work, experimental discoveries and 
the development of novel computational tools. This quest 
already gave rise to many fruitful developments, as rec- 
ognized by the 2021 Nobel Prize in Physics awarded to 
Giorgio Parisi who has been a driving force in this field. 
In addition to understanding glass more deeply, tools and 
concepts developed for glass physics are at the core of 
deep analogies between a multitude of complex systems, 
such as neural networks in biology [16, 62], algorithms for 
constrained optimisation [63], or, more recently, machine 
learning [64, 65]. 
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